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Abstract 

Background: Interstitial lung disease is a major cause of nnorbidity and mortality in systemic sclerosis (SSc), with 
insufficiently effective treatment options. Progression of pulmonary fibrosis involves expanding populations of 
fibroblasts, and the accumulation of extracellular matrix proteins. Characterisation of SSc lung fibroblast gene 
expression profiles underlying the fibrotic cell phenotype could enable a better understanding of the processes 
leading to the progressive build-up of scar tissue in the lungs. In this study we evaluate the transcriptomes of 
fibroblasts isolated from SSc lung biopsies at the time of diagnosis, compared with those from control lungs. 

Methods: We used Affymetrix oligonucleotide microarrays to compare the gene expression profile of pulmonary 
fibroblasts cultured from 8 patients with pulmonary fibrosis associated with SSc (SSc-ILD), with those from control 
lung tissue peripheral to resected cancer (n=10). Fibroblast cultures from 3 patients with idiopathic pulmonary 
fibrosis (IPF) were included as a further comparison. Genes differentially expressed were identified using two 
separate analysis programs following a set of pre-determined criteria: only genes significant in both analyses were 
considered. Microarray expression data was verified by qRT-PCR and/or western blot analysis. 

Results: A total of 843 genes were identified as differentially expressed in pulmonary fibroblasts from SSc-ILD and/ 
or IPF compared to control lung, with a large overlap in the expression profiles of both diseases. We observed 
increased expression of a TGF-|3 response signature including fibrosis associated genes and myofibroblast markers, 
with marked heterogeneity across samples. Strongly suppressed expression of interferon stimulated genes, 
including antiviral, chemokine, and MHC class 1 genes, was uniformly observed in fibrotic fibroblasts. This 
expression profile includes key regulators and mediators of the interferon response, such as STATl, and CXCLIO, and 
was also independent of disease group. 

Conclusions: This study identified a strongly suppressed interferon-stimulated gene program in fibroblasts from 
fibrotic lung. The data suggests that the repressed expression of interferon-stimulated genes may underpin critical 
aspects of the profibrotic fibroblast phenotype, identifying an area in pulmonary fibrosis that requires further 
investigation. 
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Background 

Pulmonary fibrosis, characterised by the destruction of 
lung architecture leading to organ failure, is, together 
with pulmonary hypertension (PH), the leading cause of 
death in patients with systemic sclerosis (SSc) [1]. Intersti- 
tial lung disease is more common in SSc (SSc-ILD) than 
in any other connective tissue disease, occurring in > 70% 
of patients [2], and is most frequently associated with a 
pattern of non-specific interstitial pneumonia (NSIP) [3]. 
Despite intense research efforts, the underlying mecha- 
nisms of SSc-ILD remain largely unknown [4], and there 
are currently limited therapeutic options for this serious 
complication [2]. 

While a large number of hypothesis-driven studies 
have identified potential profibrotic mediators [4,5], transla- 
tion of these into therapeutic targets has so far been largely 
disappointing [6]. The search for more effective targets in 
lung fibrosis is now being addressed using hypothesis 
generating microarray-based strategies [7,8]. The majority 
of these studies have investigated gene expression in whole 
lung tissue samples, mostly in the idiopathic setting [5]. 
Matrix metalloproteinase (MMP) 7 [9], osteopontin [10], 
Twistl [11], and MMP19 [12], are among suggested media- 
tors identified using this strategy in idiopathic pulmonary 
fibrosis (IPF), a disease characterised by a histological pat- 
tern of usual interstitial pneumonia (UIP) [13,14]. In SSc, 
most microarray studies have been performed on skin biop- 
sies and dermal fibroblasts [7]. However, a recent study 
compared whole lung tissue and fibroblasts isolated, at the 
time of transplant, from SSc-ILD lungs with a histological 
pattern of UIP, with those from IPF and idiopathic PH. The 
investigators reported gene profiles of SSc-ILD/UIP, with ei- 
ther predominant fibrosis or PH, overlapping with profiles 
of IPF and idiopathic PH, respectively [15]. 

While the initiating factors for fibrosis development 
may vary between diseases, the progressive accumulation 
of scar tissue in the lung is characterised by common 
themes, including expanding populations of activated fi- 
broblasts, and excessive accumulation of extracellular 
matrix (ECM) proteins [5]. An important strategy to iden- 
tify potential therapeutic targets, therefore, is to define fi- 
brotic fibroblast phenotypes so as to delineate underlying 
key mechanisms for fibrosis progression. 

Here we report analysis of the transcriptome of fibro- 
blasts isolated from surgical lung biopsies at the time of 
diagnosis, from patients with well defined SSc-ILD and 
the histopathological pattern of NSIP. Although the main 
aim of this study was to compare SSc-ILD/NSIP fibroblast 
gene expression profiles with those of control lung fibro- 
blasts, we also included a small number of IPF-derived 
fibroblast lines, as a separate fibrotic group. Our study 
confirms, with a robust signature in both diseases, the ab- 
errant expression of previously reported myofibroblast 
markers and fibrosis mediators, and identifies a number of 



novel, co-expressed putative disease targets. We also ob- 
served the suppression of a large gene program, the inter- 
feron stimulated genes (ISGs), reported here for the first 
time. From the known function of some of these genes 
[16], it is possible to hypothesise that this suppressed pro- 
gram underlies key fibrotic fibroblast properties, such as 
hyper-proliferation, and apoptosis resistance. This study 
therefore identifies a potential new area for investigation 
and possible intervention in pulmonary fibrosis. 

Methods 

Patients and primary lung fibroblasts 

Primary adult pulmonary fibroblasts were cultured from 
control tissue samples of unaffected lung from patients 
undergoing cancer-resection surgery (n=10), and from 
surgical lung biopsy samples of 11 patients with pul- 
monary fibrosis (SSc-ILD n=8 and IPF n=3). Independ- 
ent reviews of the clinical and histopathologic diagnoses 
were performed and conformed to established criteria 
[17,18]. All of the SSc-ILD biopsies were characterised 
by a fibrotic NSIP pattern, and the IPF biopsies by a UIP 
pattern, based on current consensus criteria for these 
histological patterns [19]. The control tissue was histo- 
logically normal. Median age (range) was 60 (52-78) in 
controls, 48 (38-69) in SSc-ILD, and 61 (44-67) in IPF. 
The gender distribution (M/F) was as follows: controls 
6/4; SSc-ILD 2/6; IPF 2/1. Four of the SSc-ILD and two 
of the IPF patients were ex-smokers. Smoking status was 
not available for all control cases. In SSc-ILD patients, 
median (range) percent predicted FVC was 72.5% (61- 
106), median FEVl was 79% (58-92) and median DLCO 
was 50% (24-58). In IPF patients, median FVC was 70% 
(64-75), median FEVl was 66% (55-79), and median 
DLCO was 50% (35-53). Patients had not been on corti- 
costeroids or other immunosuppressants prior to surgical 
biopsy, as the biopsy was performed at the time of diagno- 
sis of the ILD pattern, prior to initiation of treatment. In- 
formed written consent was given by all subjects, and 
authorisation given by the Royal Brompton Hospital Eth- 
ics Committee. Fibroblasts were obtained from the biop- 
sies by explant culture, and cell cultures maintained, as 
previously described [20,21]. Experiments were performed 
on fibroblasts at passage 2-5. Only one sample (SI) was 
used at passage 2. There was no difference in the median 
passage number between the control (median: 4.5; range: 
3-5), SSc-ILD (median: 4; range: 2-5) and IPF (median: 4; 
range 3-5) groups. 

Microarray gene expression and analysis 

At confluence, lung fibroblasts were serum-deprived for 
42 hours (media changed at 18 hours) in the presence of 
0.1% bovine serum albumin (Sigma). Total RNA was 
harvested (Trizol, Life Technologies), quantified, and the 
integrity verified by denaturing gel electrophoresis. Samples 
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with a 28S:18S ratios of approximately 2:1 were accepted 
for further analysis by the Genomics Laboratory, CSC- 
MRC, Imperial College London, Hammersmith. RNA 
samples were prepared for chip hybridisation following 
manufacturers guidelines (Affymetrix). Hybridisation of 
cRNA to Aff)^metrix human U133Av2 chips, containing 
approximately 14,500 well characterised human genes, 
signal amplification, and data collection were performed 
using an Affymetrix fluidics station and chip reader, 
according to manufacturer s protocol Array normalisa- 
tion, using the invariant set normalisation method, and 
subsequent calculation of model-based expression values, 
was performed using DNA-Chip Analyzer (dChip) [22]. A 
list of differentially expressed genes was generated in 
dChip using fold change >2, difference in means >100, 
and p<0.05. Significance analysis of microarrays (SAM) v 
4.0 [23] was also used to determine significantly differen- 
tially expressed genes with fold change >2, difference in 
means >100, delta =1, and a false discovery rate <0.01. 
Only genes identified as differentially expressed according 
to both programs were considered as different between 
groups. Microarray data has been deposited in the Gene 
Expression Omnibus database [24], accession number 
GSE40839. Although the main aim of this study was to as- 
sess global gene expression profiles in SSc-ILD compared 
to controls, for completeness we present the comparison 
between both SSc-ILD and IPF and controls separately in 
tables. dChip software was used for data visualisation and 
hierarchical average linkage clustering using Pearson s cor- 
relation [22]. 

Functional category analysis 

Functional categories enriched in the differentially 
expressed genes were identified using the functional an- 
notation and clustering tool of the Database for Annota- 
tion, Visualisation, and Integrated Discovery (DAVID) v 
6.7 [25,26]. The probability that a Gene Ontology (GO) 
biological process term [27] is overrepresented was de- 
termined by a modified Fishers exact test, comparing 
the proportion of genes in the whole genome which are 
part of that GO term, to the proportion of the differentially 
expressed genes which are part of the same GO term, and 
was expressed as an EASE score. Clusters of overrepre- 
sented GO terms were then generated based on the simi- 
larity of differentially expressed genes assigned to each 
functional GO term. Clusters were considered significantly 
overrepresented if they contained a minimum of five GO 
terms with an EASE score of < 0.01. A summary descrip- 
tion of each cluster was generated based on the constitu- 
tive GO term names of that cluster which achieved an 
EASE score <0.05 following Benjamini-Hotchberg correc- 
tion of multiple comparisons. Only clusters with enrich- 
ment scores >3 (minus log transformed geometric mean of 
the EASE scores of the constitutive terms, equivalent to 



average EASE score=0.001) were selected. The open access 
database INTERFEROME [28] was used to identify differ- 
entially expressed genes which have been shown experi- 
mentally to be regulated by interferons. 

qRT-PCR 

RNA was extracted using the RNeasy® Mini kit (Qiagen) 
according to manufacturers instructions. Samples were 
quantified and quality assessed using the NanoDrop 
spectrophotometry system (Thermo Scientific). Com- 
plementary DNA (cDNA) was synthesised from 500 ng 
RNA in a 20 [i\ reaction using the QuanTitect® reverse 
transcription kit (Qiagen). Expression levels were mea- 
sured using a Rotor Gene 6000 (Corbett) in 10 [A reac- 
tions containing 2 [A cDNA (five-fold dilution), 1 x 
SensiMix™ SYBR NO-ROX (BioUne), and 0.5 [iM of each 
forward and reverse primer (Table 1). PCR conditions 
were: 10 minutes at 95°C, followed by 40 cycles of 
10 seconds at 95°C, 15 seconds at 57°C, and 5 seconds 
at 72°C. All reactions were performed in duplicate and 
non-template controls were included for each gene. 
Standard curves were generated for each gene studied 
using seven two-fold serial dilutions, high standard of 
1x10^ copies/(il, of primer set amplicons generated from 
cDNA. Threshold cycle was manually determined at a 
fixed value of 10'°'^ and the template quantity calculated 
using Rotor Gene 6000 Series Software 1.7 (Corbett). 
Expression levels were normalised to YWHAZ and 
HPRTl [29,30]. 

Western blot analysis 

Following pre-incubation for 24 hours in serum-free 
media (DMEM, 0.1% BSA, penicillin/streptomycin), pul- 
monary fibroblasts from healthy controls, SSc-ILD, and 



Table 1 qRT-PCR primers 



Gene 


Forward primer 


Reverse primer 


Normalisation genes 


HPRTl 


TGACACTGGCAAAACAATGCA 


GGTCC^CACCAGCAAGCT 


YWHAZ 


AC^GGTACATOTGGOTCAA 


CCGCCAGGACAAACCAGTAT 


Genes of interest 




CXCLIO 


GAAAGCAGTOGCAAGGAAAG 


ATCOTGGAAGCACTGCATC 


IDl 


CCAGAACCGCAAGGTGAG 


GGTCCCTGATGTAGTCGATGA 


IFITMl 


TOTOAACTGGTGCTGTCT 


ATGAGGATGCCCAGAATCAG 


ILll 


CCTGTGGGGACATGAACTGT 


AGGGTCTGGGGAAACTCG 


IRFl 


CAGCCCAAGAAAGGTCCTC 


TOAACGGTACAGACAGAGCA 


N0X4 


CTGCTGACGTOCATGmC 


CGGGAGGGTGGGTATCTAA 


Serpine 
1 


GGAAAGGCAACATGACCAG 


CAGGTOTCTAGGGGOTCC 


STATl 


GGATCAGCTGCAGAACTGGT 


mCTGTCCAATOCTCCAA 



Shown are the primers, written 5' to 3', used for measuring expression levels 
by qRT-PCR for validation of the microarray results. All primer pairs except for 
YWHAZ, IRFl, and IDl, are intron spanning. 
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IPF (n=3 for each phenotype), were cultured for a fur- 
ther 24 hours in fresh serum-free media. Cells were 
lysed and western blot analysis was performed using the 
following primary antibodies: CTGF and STATl (Santa 
Cruz Biotechnology); aSMA (Dako); IFITMl-3, ISG15 
and GAPDH (Abeam); IRF-1 (Cell Signaling Technol- 
ogy); horseradish peroxidise conjugated secondary anti- 
bodies (Dako and Cell Signaling Technology); and ECL 
detection (Amersham). 

Results 

Gene expression profiles of fibrotic lung fibroblasts: 
approximately two-thirds of differentially expressed 
genes are down-regulated 

Using an Affymetrix platform (U133Av2), we determined 
basal (serum free) global gene expression levels in fibro- 
blasts prepared from lung tissue of 8 patients with SSc- 
ILD and 10 control lungs. As a further comparison we 
also included 3 fibroblast cultures from lung tissue of 
IPF patients. Unsupervised hierarchical cluster analysis 
of samples and genes resulted in an overall separation of 
fibrotic samples from controls (Figure 1). Two of the 
SSc-ILD samples clustered among the normal controls, 
demonstrating recognised fibrotic fibroblast sample 
heterogeneity. 

When stringent criteria were applied, as described in 
the methods section, 478 and 744 probe sets (probes in 
future), equivalent to 360 and 547 genes, displayed dif- 
ferential expression by at least two-fold in SSc-ILD and 
IPF fibroblasts, respectively, compared to control fibro- 
blasts. In SSc-ILD fibroblasts, 125 probes (99 genes) 
showed significantly higher, and 353 probes (261 genes) 
lower expression levels compared with controls (Additional 
file 1). In IPF fibroblasts, 239 and 505 probes (181 and 366 
genes) had significantly higher and lower expression, re- 
spectively, compared with control fibroblasts (Additional 
file 2). Thus, approximately two thirds of differentially 
expressed genes were expressed at lower levels in fibro- 
blasts from the two disease groups compared with con- 
trols. The sets of underexpressed genes in the disease 
groups, compared to normal controls, also contained the 
most significantly differentially expressed genes. Table 2 
lists the 20 most significantly differentially expressed, and 
highest ranking (fold change) genes in SSc-ILD and IPF 



samples, separately, versus controls. It should be noted 
here that, while 379 out of a total of 843 probes were 
differentially expressed in both SSc-ILD and IPF vs. 
controls, most genes falling outside this overlap were 
very nearly significant in the other disease group, i.e. 
followed similar trends without meeting cut-off criteria. 
When the two disease groups were compared directly, 
only 8 probes (7 genes) were significantly differentially 
expressed (Additional file 3). The 843 probes differen- 
tially expressed in SSc-ILD and/or IPF, which also in- 
clude the genes differentially expressed between these 
two groups, are listed in Additional file 4. Within the 
pooled fibrotic samples (IPF and SSc-ILD), no significant 
difference was observed in the 843 probes according to 
smoking status (SAM analysis, data not shown). 

Expression levels of a subset of genes identified by the 
present microarray analysis were verified by qRT-PCR in 
the microarray RNA samples demonstrating good cor- 
relation between these two techniques (Figure 2). Pro- 
tein levels for six differentially expressed genes, fibrosis 
related genes; connective tissue growth factor (CTGFJ, 
and alpha-smooth muscle actin (aSMA), and interferon 
stimulated genes (ISG); signal transducer and activator 
of transcription 1 (STATl), IFITMl-3 (all three isoforms 
are detected by this antibody), ISG 15 and IRF-1, were 
determined by western blot analysis in independent 
preparations of additional SSc-ILD, IPF, and control fi- 
broblasts (Figure 3). 

Functional groups over-represented among differentially 
expressed genes 

As predicted from the similarities between the two disease 
groups with regards to differentially expressed genes, gene 
ontology analysis revealed a major overlap in overrepre- 
sented functional groups. Enriched functional groups 
representing six broad categories (GO clusters) were iden- 
tified among the genes with higher expression in disease fi- 
broblasts compared with normal controls (see Table 3 for a 
summary of enriched clusters of GO terms, and Additional 
file 5 for full GO analysis): anatomical structure develop- 
ment, regulation of cell cycle, response to stress and 
wounding, regulation of apoptosis, cell migration, and 
smooth muscle contraction. Among underexpressed/ 
downregulated genes, GO clusters included: inflammatory 
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Figure 1 Unsupervised clustering of samples based on full microarray probe set. The sample dendrogram resulting from hierarchical 
clustering using all 22 K probes, shows clustering of samples by phenotype: control (C, green bar), SSc-ILD (S, orange bar), IPF (U, red bar). 
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Table 2 Most differentially expressed genes in SSc-ILD and IPF compared to controls 



Genes overexpressed in SSc-ILD vs. Control 


Probe set ID 


Control mean 


SSc-ILD mean 


Fold change 


p value 


Inhibitor of DNA binding 1, dominant negative lielix-loop-lielix protein 


208937_ 


_s_at 


25.5 


917.5 


36.1 


0.00078 


lnterleul<in 1 1 


206924_ 


_at 


23.6 


717.8 


30.4 


0.015 


Inliibitor of DNA binding 3, dominant negative lielix-loop-lielix protein 


207826. 


_s_at 


27.7 


603.2 


21.8 


0.00051 


Tetraspanin 13 


217979. 


_at 


37.9 


533.9 


14.1 


0.0033 


Elastin 


212670. 


_at 


43.7 


396.9 


9.1 


0.0021 


Xylosyltransferase 1 


213725. 


_x_at 


29.0 


255.4 


8.8 


0.0024 


Serpin peptidase inliibitor, clade E, member 1 


202628. 


_s_at 


329.8 


2473.2 


7.5 


0.0022 


Basic lielix-loop-lielix family, member e40 


201170. 


_s_at 


44.3 


256.1 


5.8 


0.0014 


Connective tissue growth factor 


209101. 


_at 


467.9 


2637.1 


5.6 


0.00068 


Solute carrier family 7, member 5 


201195. 


_s_at 


56.1 


294.2 


5.3 


0.00022 


Tropomyosin 1 (alpha) 


206116. 


_s_at 


317.5 


1619.6 


5.1 


0.0032 


Phosphoribosyl pyrophosphate synthetase 1 




_b_d L 


59.0 


283.4 


4.8 


0.0036 


Inhibin, beta A 


1 1 r\Q 1 1 
z 1 U J 1 1 . 


_s_at 


1 /I Q Q 


(^QQ 1 
Doo. 1 


A Q 
4.0 


U.UU 1 z 


Growth arrest and DNA-damage-inducible, beta 


207574. 


_s_at 


68.1 


306.6 


4.5 


0.0023 


Coiled-coil domain containing 99 


221685. 


_s_at 


85.1 


373.3 


4.4 


0.0018 


Cadherin 2, type 1, N-cadherin (neuronal) 


203440. 


_at 


105.3 


433.0 


4.1 


0.00095 


Desmoplakin 


200606. 


_at 


80.5 


306.1 


3.8 


0.0016 


Insulin-like growth factor binding protein 3 


212143. 


_s_at 


408.6 


1489.1 


3.6 


0.0068 


Microtubule associated monoxygenase, calponin and LIM domain 
containing 2 


212473. 


_s_at 


184.7 


667.5 


3.6 


0.0039 


Prostaglandin-endoperoxide synthase 1 


215813. 


_s_at 


102.4 


362.0 


3.5 


0.0083 


Genes underexpressed in SSc-ILD vs. Control 


Probe set ID 


Control mean 


SSc-ILD mean 


Fold change 


p value 


Chemokine (C-X-C motif) ligand 10 


204533. 


_at 


771.2 


19.2 


-40.1 


0.00034 


Chemokine (C-X-C motif) ligand 1 1 


210163. 


_at 


179.9 


5.0 


-36.0 


0.0027 


Flavin containing monooxygenase 2 (non-functional) 


211726. 


_s_at 


530.4 


15.7 


-33.7 


0.017 


Interferon-induced protein with tetratricopeptide repeats 2 


217502. 


_at 


707.1 


26.1 


-27.1 


0.0096 


Vascular cell adhesion molecule 1 


203868. 


_s_at 


835.2 


32.2 


-26.0 


0.0049 


Bone marrow stromal cell antigen 2 


201641. 


_at 


315.8 


12.5 


-25.3 


0.0081 


Radical S-adenosyl methionine domain containing 2 


213797. 


_at 


333.8 


13.3 


-25.1 


0.0039 


Interferon-induced protein 44-like 


204439. 


_at 


370.7 


15.4 


-24.1 


0.00033 


Interferon-induced protein with tetratricopeptide repeats 1 


203153. 


_at 


1 744.2 


82.6 


-21.1 


0.000039 


2',5'-oligoadenylate synthetase 1, 40/46 kDa 


205552. 


_s_at 


374.6 


18.3 


-20.5 


0.0014 


Complement factor B 


202357. 


_s_at 


837.0 


42.6 


-19.6 


0.0036 


Interferon-induced protein with tetratricopeptide repeats 3 


204747. 


_at 


985.0 


61.9 


-15.9 


0.00031 


Chromosome 10 open reading frame 10 


209183. 


_s_at 


208.8 


13.6 


-15.3 


0.0062 


Myxovirus resistance 1, interferon-inducible protein p78 (mouse) 


202086. 


_at 


1361.9 


91.2 


-14.9 


0.00012 


Receptor (chemosensory) transporter protein 4 


219684. 


_at 


196.4 


13.3 


-14.8 


0.000091 


Chemokine (C-C motif) ligand 1 1 


210133. 


_at 


529.8 


36.2 


-14.6 


0.0021 


Retinoic acid receptor responder (tazarotene induced) 3 


204070. 


_at 


239.4 


17.2 


-13.9 


0.00024 


Alcohol dehydrogenase IB (class 1), beta polypeptide 


209613. 


_s_at 


268.9 


19.8 


-13.6 


0.011 


Secreted and transmembrane 1 


213716. 


_s_at 


285.0 


22.0 


-13.0 


0.0012 


Interferon, alpha-inducible protein 6 


204415. 


_at 


1196.6 


93.7 


-12.8 


0.00043 


Genes overexpressed in IPF vs. Control 


Probe set ID 


Control mean 


IPF mean 


Fold change 


p value 


Interleukin 1 1 


206924. 


_at 


23.6 


2374.9 


100.6 


0.0019 


Inhibitor of DNA binding 1, dominant negative helix-loop-helix protein 


208937. 


_s_at 


25.5 


752.8 


29.6 


0.011 
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Table 2 Most differentially expressed genes in SSc-ILD and IPF compared to controls (Continued) 



Tetraspanin 13 
NADPH oxidase 4 

Inhibitor of DNA binding 3, dominant negative lielix-loop-lielix protein 

Pliospliolamban 

Elastin 

Xylosyltransferase I 
Galanin prepropeptide 
Cytol<ine receptor-! il<e factor 1 
Calponin 1, basic, smootli muscle 
Follistatin-like 3 
CTP synthase 

Endothelial cell-specific molecule 1 
Cadherin 6, type 2, K-cadherin (fetal kidney) 
Proenkephalin 

Adhesion molecule with Ig-like domain 2 
NUAK family, SNFl-like kinase, 1 
Tropomyosin 1 (alpha) 
Inhibin, beta A 

Genes underexpressed in IPF vs. Control 



and immune response, response to biotic stimulus, regula- 
tion of apoptosis, regulation of cell migration, regulation of 
cell proliferation, and regulation of I-kB/NF-kB cascade. 
By far the most enriched functional groups, compared 
with control cells, in both SSc-ILD and IPF lung 



.at 


37.9 


1 039.1 


27.4 


0.0052 


.at 


1 2.3 


323.6 


26.4 


0.016 


.s_at 


27.7 


603.4 


21 .8 


0.0025 


.s_at 


25.7 


460.1 


1 7.9 


0.035 


.at 


43.7 


766.9 


1 7.6 


0.0026 


.x_at 


29.0 


443.6 


1 5.3 


0.034 


.at 


1 8.4 


254.4 


1 3.8 


0.014 


.at 


23.4 


31 9.3 


1 3.6 


r\ r\r\no 

u.uuyo 


.at 


83.0 


1 048.0 


1 2.6 


0.0024 


.s_at 


32.5 


404.8 


1 2.4 


n, r\r\/\ o 

(J.(J(J4o 


at 


39.5 


454.7 
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fibroblasts, are genes involved in immune system pro- 
cesses and in microbial/viral defence, which are strongly 
suppressed in both disease groups. These genes are also 
among the most significantly differentially expressed 
genes in this study. 
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Shown are the top 20 genes, based on fold change, over and under expressed in SSc-ILD and IPF compared to controls. Where more than one probe set 
corresponding to the same gene were present in the top 20, the probe with the greatest fold change or most significant p-value, is shown. 
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■ Microarray ^ qRT-PCR 

Figure 2 Quantitative RT-PCR confirmation of microarray results. Expression levels of eight genes selected from the microarray data was 
measured by qRT-PCR in thirteen of the samples used in the microarray. For each sample the microarray data is plotted on the left-hand axis, and 
the qRT-PCR results plotted on the right-hand axis. qRT-PCR expression levels were normalised to YWHAZ and HPRTl. 



Cluster analysis of differentially expressed genes - ISGs; 
a major group suppressed in fibrotic lung fibroblasts 

To visualise the differential expression across samples 
and to identify co-expressed genes, average linkage clus- 
ter analysis was performed using expression data for the 
843 probes sets which displayed differential expression 
in at least one of the two comparisons: SSc-ILD vs. con- 
trol, and IPF vs. control (Figure 4A). Within this set are 



also the eight probes differentially expressed when SSc- 
ILD and IPF samples were compared directly. Parts of 
identified gene clusters were selected to illustrate co- 
expression among upregulated (Figure 4, Panels B-F) 
and downregulated (Panels H-J) genes, and also to high- 
light different patterns of sample heterogeneity. Overall, 
a heterogeneous expression pattern was observed in the 
upregulated genes, whereas downregulated genes had 
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Figure 3 Western blot confirmation of microarray results. Protein expression levels of six significantly differently expressed genes selected 
from the microarray data were visualised by western blot in fibroblast samples independent to those used in the microarray (n=3 for each 
group). A) Protein expression of fibrosis related genes, CTGF and aSMA; and key regulator of interferon response, STATl. B) Protein expression of 
interferon stimulated genes (ISG): IFITMl-3, ISG15 and IRF-1. 



more uniform expression patterns across samples for the 
majority of genes in both disease groups. Among 
upregulated genes, a TGF-p response signature [31,32] 
including genes encoding for fibrosis mediators and 
myofibroblast markers, such as SERPINEl {PAH), con- 
nective tissue growth factor, smoothelin, and transgelin 
(SM22) is prominent (Panel B). Co-expressed with these 
are strongly upregulated genes: growth arrest and DNA- 
damage inducible |3 {GADD45), xylosyltransferase 1 
(XYLTl), N-cadherin, and elastin, with potential roles in 
the fibrotic disease process. Groups B, C and D all con- 
tain genes involved in contraction and migration, how- 
ever, the degree of heterogeneity between samples differ 
between these groups: in group B, the majority of fi- 
brotic samples have elevated expression of smoothelin 
and transgelin compared with controls; in group C, 
fewer fibrotic samples, 7 out of 11, have enhanced levels 
of a2 smooth muscle actin {ACTA2) expression; and in 
group D, all three IPF samples, but only 3 out of 8 SSc- 
ILD samples, have elevated expression of calponin 1 and 
actin gamma 2 smooth muscle {ACTG2), This may indi- 
cate different degrees of contractile/migratory pheno- 
types among these fibrotic cell preparations. Panel E 
contains IDl and IDS, which are in the top 20 differen- 
tially regulated genes in both disease groups, and are 
upregulated in most of the fibrotic samples. Group F de- 
picts a cluster of co-expressed cell-cycle associated genes, 
including cyclins and T0P02, which exhibit heterogeneous 
expression in both disease groups. Panel G illustrates an 
area with less clustering, which however includes possible 
disease specific genes, e.g. Secreted protein, acidic, 
cysteine-rich {SPARC) (IPF, Table 2B) and desmoplakin 
(SSc-ILD) (Additional file 4). Desmoplakin is among the 
top 20 most upregulated genes in SSc-ILD with an elevated 
expression in 7 out of the 8 SSc-ILD fibroblast lines, but 
with low expression in the three IPF cell preparations and 
in all controls. Desmoplakin is part of the desmosome 
complex which forms tight cell-cell contacts [33], and its 



enhanced expression in SSc-ILD fibroblasts may define a 
different pathogenesis and cell origin. Panels H-J shows the 
marked suppressed expression of interferon stimulated 
genes (ISGs), such as antiviral genes (Group H), chemokine 
(Group I) and MHC class I genes (Group J). This cluster 
also includes key regulators of the interferon gene pro- 
gram, STATl, interferon regulatory factor 1 {IRFl) and 
interferon regulatory factor 7 {IRF7), as well as chemokine 
(C-X-C motif) ligand 10 (CXCLIO/IPIO), one of the most 
strongly suppressed genes in the study, and the most 
strongly repressed chemokine in both disease groups 
(Panel H and I). 

Comparing the SSc-ILD and IFP fibroblast gene 
expression profiles with the interferome 

Since we observe a clear fibrosis/TGF-|3 signature to- 
gether with a strongly suppressed ISG program, and 
there is a well-documented antagonistic relationship be- 
tween TGF-|3 and interferon signalling in the fibrosis lit- 
erature [34], we next interrogated the Interferome, a 
database of interferon regulated genes (IRGs) reported 
in the literature [28]. This database includes 1996 hu- 
man IRGs, of which 1581 are induced and 415 repressed 
by interferons. In our study, out of the 99 and 181 
overexpressed genes in SSc-ILD and IPF fibroblasts, re- 
spectively, 40 (40.4%) and 55 (30.1%) genes were in the 
Interferome database. Among our underexpressed genes, 
out of a total of 261 and 366 in SSc-ILD and IPF fibro- 
blasts, respectively, 134 (51.3%) and 173 (74.3%) genes 
were in the Interferome database. The genes overlapping 
with IRGs in this database are listed in Additional file 6. 
The comparison revealed that many of the TGF-p re- 
sponsive genes upregulated in our microarray data set 
are indeed genes repressed by IFNs. 

Discussion 

Previous studies have shown that fibroblasts isolated 
from SSc-ILD [35] and IPF [5] lungs, while displaying 
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Table 3 Over-represented functional terms 


GO cluster 


Description 


Enrichment 
score 


Genes overexpressed in SSc-ILD vs. Control 


1 


Anatomical structure development 


5.50 


2 


Regulation of cell cycle 


4.23 


3 


Response to stress and wounding 


3.33 


Genes underexpressed In SSc-ILD vs. Control 




1 


Inflammatory response 


11.53 


2 


Regulation of cell proliferation 


5.16 


3 


Regulation of biological process 


5.12 


4 


Inflammatory response/chemotaxis 


4.81 


5 


Regulation of cell migration 


4.33 


6 


Response to external stimulus 


4.19 


7 


Regulation of apoptosis 


3.78 


8 


Inflammatory and immune response 


3.77 


9 


Response to biotic stimulus/ion homeostasis 




10 


Regulation of I-kB kinase/NF-KB cascade 


D.D/ 


Genes overexpressed in IPF vs. Control 




1 


Anatomical structure development/ 
neurogenesis 


4 94 


2 


Regulation of apoptosis 


D.yJ 


3 


Cell migration/neurogenesis 


3.81 


4 


Regulation of cell motion 


3.45 


5 


Response to wounding/tissue development 


3.17 


6 


Smooth muscle contraction/Blood 
circulation 


3.03 


Genes underexpressed in IPF vs. Control 




1 


Immune response 


12.65 


2 


Response to virus, bacteria and LPS 


8.66 


3 


Positive regulation of biological process and 
cell death 


5.29 


4 


Negative regulation of biological process 
and cell death 


4.91 


5 


Regulation of immune system and 
developmental process 


4.08 


6 


Inflammatory response/chemotaxis 


3.98 


7 


Regulation of cell migration and adhesion 


3.87 


8 


Inflammatory and humoral immune 
response 


3.63 


9 


Anatomical structure development 


3.39 


10 


Response to stimulus and I-kB kinase/NF-KB 
cascade 


3.19 


11 


Response to extracellular stimulus and 
oxidative stress 


3.18 



Using the DAVID functional annotation tool, genes over and under expressed 
in SSc-ILD and IPF compared to controls were clustered according to Gene 
Ontology (GO) biological process terms. Shown are the summary descriptions 
and enrichment scores of the sets of enriched GO terms within each GO 
cluster with an enrichment score >3. 



substantial heterogeneity, are generally more proliferative, 
migratory, resistant to apoptosis, and ECM producing, 
than control lung fibroblasts [36]. A greater proportion 
display elevated aSMA expression, also enhancing their 
contractility. These features are all consistent with the 
fibroblast accumulation and scar tissue formation observed 
in fibrotic lung. The so called myofibroblast phenotype is 
maintained over several passages in culture [37], suggesting 
an underlying epigenetic regulation, whether established 
and maintained in local cells within the chronic disease set- 
ting, or supporting the phenotype of a specialised infiltrat- 
ing wound healing cell type [38]. Regardless of cell origin, 
this feature enables in vitro studies of mechanisms under- 
lying the fibrotic fibroblast phenotype. 

In this study we observed a high number of differen- 
tially expressed genes between SSc-ILD and/or IPF de- 
rived fibroblasts, compared with controls. There was a 
large overlap between expression profiles of SSc-ILD and 
IPF fibroblasts, suggesting several common pathways at 
this stage of the two diseases. Indeed, a direct comparison 
demonstrated only seven genes with significantly different 
expression between the two disease groups. Caution 
should be applied, however, when interpreting data from 
this direct comparison since only three IPF samples were 
included here. It is possible that more genes with differen- 
tial expression between the two disease groups would be 
identified with a larger numbers of IPF fibroblasts sam- 
ples. Such studies would also be required to verif)^ the 
sporadic observations made in this study, such as the 
elevated expression levels of desmoplakin in SSc-ILD, but 
not in IPF fibroblasts. Therefore, we stress that the main 
objective of the study presented here was to gain an over- 
view of potential SSc-ILD target genes, and the study was 
not designed to detect differences between different fi- 
brotic entities. 

Among significantly upregulated genes in both SSc-ILD 
and IPF/UIP fibroblasts, we identify a recognised fibrosis 
signature, including smooth muscle actin {ACTA2), CTGF, 
and PAIl {SERPINEl), along with genes more recently as- 
sociated with lung fibrogenesis, including IDl, IDS, ILll, 
and N0X4. Inhibitor of DNA binding 1 and 3 [IDl and 
ID3)y are target genes of bone morphogenetic proteins, and 
control cell differentiation by dominant negative inhibition 
of helix-loop-helix transcription factors [39]. Chambers 
et al. described upregulation of IDl in lung fibroblasts in 
response to TGF-|3, and linkage to the smooth muscle 
phenotypic switch [31]. Increased IDl and IDS gene ex- 
pression was also observed in lung tissue and fibroblasts 
from patients with SSc-ILD by Hsu and colleagues, [15]. 
The ID proteins are overexpressed in many cancers, 
controlling cell growth and apoptosis, and have been sug- 
gested as a therapeutic target [39,40]. In fibroproliferative 
stages of fibrosis, high ID1/ID3 expression could maintain 
fibroblasts in a dedifferentiated, hyperproliferative, and 
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Figure 4 Gene expression profiles of control, SSc-ILD, and IPF pulmonary fibroblasts. A) Supervised hierarchical clustering of probes. The 
gene set used comprised the 843 probes with a fold change >2, difference in means >100, and p<0.05 (dChip) and FDR<0.01 (SAM), in both the 
SSc-ILD or IPF samples, compared to controls, which includes also 8 probes differentially expressed between SSc-ILD and IPF samples when 
compared directly Each column corresponds to an individual sample (C= control, S=SSc-ILD, U=IPF), and each row corresponds to an individual 
probe. The coloured bars to the right of the heatmap identify the location of the insets displayed in B-J. Panels B-J were selected to highlight 
sample heterogeneity and/or genes co-expressed with high scoring differentially expressed genes. 



apoptosis resistant state. Another significantly upregulated 
gene, though somewhat variably in SSc-ILD fibroblasts, is 
N0X4, encoding a member of the NADPH oxidase (NOX) 
proteins which generate superoxide by electron transfer to 
oxygen [41]. Through its involvement in TGF-|3-induced 
fibroblast differentiation into myofibroblasts, two recent 
studies have suggested a role for NOX4 in IPF, where it is 
found to be overexpressed [41]. Specific NOX4 inhibitors 
are now being developed as possible antifibrotic agents 
[42]. Another potential novel target gene, GADD45B, 
highly induced and coexpressed with the fibrosis related 
genes in this study, is a pro- survival factor associated 
with stress-resistant tumours [43], and has been found 
to be upregulated in SSc skin biopsies [44]. 

The most striking of our observations is the strongly 
and uniformly repressed ISG profile in fibrotic fibro- 
blasts. Among these are genes coding for: antiviral 



protein myxovirus (influenza virus) resistance 1 (MXl), 
interferon gamma inducible protein pl6 {IFI16), XS'- 
oligoadenylate synthetase 1 {OASl)^ the chemokines 
chemokine (C-X-C motif) ligand 10 (CXCLIO/IPIO) and 
chemokine (C-C motif) ligand 11 {CCLll), and antigen 
presenting MHC I molecules. Key transcriptional regula- 
tors of this program, including IRFl, IRF7, and STATl 
are also suppressed. To our knowledge, the reduction in 
expression of a large set of immune response/interferon 
related genes has not previously been described in either 
SSc-ILD or IPF derived fibroblasts, or other fibrotic lung 
diseases. Global gene expression was recently evaluated 
in scleroderma whole skin biopsies and matched dermal 
fibroblasts. While 26 genes were differentially expressed 
in both scleroderma whole skin and fibroblasts com- 
pared to controls, nine were found to be discordant [45]. 
Interestingly, the majority of the discordant genes were 
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upregulated in skin biopsies but downregulated in fibro- 
blasts, including the ISG genes we identify as repressed 
in lung fibroblasts, such as MXl, IFI16, intercellular ad- 
hesion molecule 1 {ICAMl), and OASl, While the au- 
thors argue that the discordant genes are an indication 
of a scarce representativeness of skin fibroblast gene ex- 
pression in vivo [45], in light of our data, it is likely that 
they were in fact observing the same phenomenon 
reported in the current study. As these are mainly inter- 
feron regulated genes, the discordance between fibroblasts 
and whole tissue gene expression could be explained by 
the known increase of immune cell populations in SSc 
skin, likely to be overexpressing immunoregulatory genes. 
The upregulation of MXl, IRF7 and STATl in PBMC 
from SSc patients would support this notion [46,47]. The 
downregulation of the ISG program both in SSc-ILD and 
in IPF fibroblasts suggests that this phenomenon relates 
to a local fibroblast specific, rather than systemic, pro- 
fibrotic process, perhaps underpinned by a general suscep- 
tibility for tissue fibrosis, common to both diseases. 

Additional support for our finding of an aberrantly 
regulated ISG program in fibrotic lung fibroblasts, comes 
from work on two of the signature genes from this group, 
CXCLIO (IPIO) and STATL IPIO levels were found to be 
downregulated in IPF lung fibroblasts by Keane et al. [48] . 
More recently. Coward et al. have shown, again in IPF fi- 
broblasts, that epigenetic dysregulation involving both his- 
tone deacetylation and hypermethylation is responsible for 
targeted repression of IPIO [49]. By contrast, Hsu et al. 
did not report a significant difference in expression in iso- 
lated lung fibroblasts [15], a difference which may relate 
to experimental design, as discussed below. Further sup- 
port of a role of suppressed ISGs in pulmonary fibrosis 
comes from two animal models. Mice deficient in IPIO 
[50], and in STATl [51], displayed enhanced susceptibil- 
ity to pulmonary fibrosis. Based on these studies, the 
suppressed ISG program indentified in fibroblasts iso- 
lated from SSc-ILD and IPF lung, as presented here, 
would support enhanced lung fibrosis progression 
through promoting fibroblast proliferation, migration, 
and apoptosis resistance. Interestingly, IFN-y treatment 
in IPF has failed to show a benefit [52], and there was a 
suggestion of worse pulmonary outcomes in a study in- 
vestigating treatment of SSc patients with IFN-a com- 
pared with placebo [53], with the latter providing 
indirect support for IFN-related mechanisms involved 
in organ-specific SSc compUcations. There could be 
several possible explanations for these disappointing re- 
sults, including unexpected adverse effects through 
circulating cell populations. Activation of pathways 
downstream of systemically administered interferons is 
likely to have different direct and indirect effects de- 
pending on the cell type and tissue location. The find- 
ings shown here add important information to this 



complexity, and need to be investigated in future de- 
tailed mechanistic cell and animal studies. 

The signatures observed in our study are remarkably 
strong, both in terms of fold difference and statistical 
significance. One possible reason for this is that lung tis- 
sue samples were obtained from biopsies at the time of 
diagnosis, when the disease may be at a relatively early 
or active wound healing stage, and when perhaps fibro- 
blast proliferation (accumulation) and elastin synthesis, 
rather than contraction and collagen remodelling, dom- 
inate. This is in contrast to the study by Hsu et al, in 
which gene expression profiles were investigated in SSc 
lung tissue and fibroblasts from transplant, and therefore 
possible end-stage, material, also noted by the authors as 
a potential limitation [15]. Another possibility relates to 
differences in in vitro culture condition as we employed 
serum free media before harvest, similarly to others, in- 
cluding Coward et al. [49], as opposed to in low serum 
(0.5%) as applied by e.g. Hsu et al. [15]. Serum with- 
drawal, a form of cellular stress, may evoke the clear 
differential expression profiles observed in our study. 
Whereas in the fibrotic fibroblasts an anti-apoptotic sur- 
vival gene program is maintained, which may be a result 
of the suppressed ISGs, this gene program may not be 
subject to repression in the normal fibroblasts. 

While an accepted source of control tissue in studies 
of ILD [10,12], the use of control fibroblasts from cancer 
resected specimens, rather than from healthy control 
subjects (not available for this study), represents a po- 
tential limitation. Although obtained from areas of lung 
with normal histological appearance, differential gene 
expression in fibroblasts derived from lungs in which 
cancer has developed cannot be excluded. However, sig- 
nificant gene expression differences were observed in 
lung cancer associated fibroblasts compared to matched 
fibroblasts from areas of normal lung from the same 
patient, suggesting that the cancer associated phenotype 
of lung fibroblasts is regionally limited to the cancer 
stroma [54]. A further possible source of bias in gene ex- 
pression is smoking history. Smoking has been shown to 
be associated with interstitial fibrotic changes [55], and 
is itself likely to cause changes in the expression of 
certain genes. We observed no significant differences 
among differentially expressed genes according to smok- 
ing status in the pooled fibrotic samples, suggesting that 
the observed changes were related to the fibrotic lung 
disease itself rather than to smoking. However, as sub- 
group numbers were small, and it was not possible to 
separately analyse the two fibrotic lung diseases, further 
studies are needed to carefully assess the contribution of 
smoking to gene expression changes in the context of fi- 
brotic lung disease. 

While the general hypothesis, that a repressed inter- 
feron stimulated gene program at least in part underpins 
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the fibrotic fibroblast phenotype, will be tested in future 
studies, it is interesting to note the same phenomenon 
in several other clinical settings where hyperplasia and 
apoptosis resistance are key features; certain viruses, in- 
cluding high-risk human papillomaviruses (HPV), have 
evolved a mechanism to down-regulate ISGs in host 
cells as an immunoevasive strategy [56], and persistent 
HPV infection may lead to cervical cancer development; 
breast cancer metastasis is promoted by IRF7 silencing 
[57]; fibroblasts from patients with Li-Fraumeni syn- 
drome become spontaneously immortalised through the 
downregulation of interferon pathway genes [58]. Con- 
versely, IRFl expression reverts the phenotype of 
oncogenically transformed fibroblasts [59], and IRF-1 
enhancing drugs with tumour suppressing properties 
are currently being developed [60]. Many similar exam- 
ples in the literature lead to questions about whether fi- 
brosis is a pre-cancerous state [61]. The repressed, or 
aberrantly regulated, fibroblast specific interferon re- 
sponse network may therefore be a common necessary 
determinant allowing lung fibrosis progression to occur. 

In summary, in this study comparing gene expression 
profiles of fibroblasts explanted from fibrotic lung tissue 
(SSc-ILD and IFF), with control fibroblasts from areas of 
normal lung, we observe: an overall elevated expression 
of previously reported fibrosis associated genes, with 
marked heterogeneity across samples; differentially regu- 
lated myofibroblast markers which correlate with the ex- 
pression heterogeneity between samples; and a strongly 
suppressed interferon stimulated gene program, uni- 
formly present across fibrotic samples. This suppressed 
gene program displays both the greatest significance and 
largest fold differences in expression in our data set. Simi- 
larly to functional findings in parallel fields, particularly 
cancer, this group of genes, and the suppression of their 
expression, could explain essential aspects of the 
profibrotic fibroblast phenotype. This hypothesis will need 
to be tested by future studies, with particular focus on epi- 
genetic silencing as a potential underlying mechanism. 



Additional files 



Additional file 4: Summary of differentially expressed genes in SSc- 
ILD and IPF. Word file, .txt extension. This data set contains all of the 
genes up- or down- regulated in fibroblasts from at least one disease 
group compared to control fibroblasts. Included are p-values from dChip 
analysis and q-values from SAM analysis. 

Additional file 5: Functional annotation clustering analysis. Excel 
file, .xlsx extension. Using the DAVID functional annotation tool, genes 
over and under expressed in SSc-ILD and IPF compared to controls were 
clustered according to Gene Ontology biological process terms. Shown 
are the enriched terms within each annotation cluster with an EASE score 
threshold of < 0.01, and an initial group membership of 5. 

Additional file 6: Differentially expressed genes present in 
Interferome. Excel file, .xlsx extension. The data set contains all of the 
genes differentially expressed genes in SSc-ILD and IPF fibroblasts 
compared to control fibroblasts which have been shown experimentally 
to be regulated by interferons listed in the INTERFEROME database. 
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Additional file 1: Genes differentially expressed in SSc-ILD. Word 
file, .txt extension. This data set contains all of the genes up- or down- 
regulated in SSc-ILD fibroblasts compared to control fibroblasts. Included 
are p-values from dChip analysis and q-values from SAM analysis. 

Additional file 2: Genes differentially expressed in IPF. Word file, txt 
extension. This data set contains all of the genes up- or down- regulated 
in IPF fibroblasts compared to control fibroblasts. Included are p-values 
from dChip analysis and q-values from SAM analysis. 

Additional file 3: Genes differentially expressed between IPF and 
SSc-ILD. Word file, .txt extension. This data set contains all of the genes 
up- or down- regulated in IPF fibroblasts compared to SSc-ILD fibroblasts. 
Included are p-values from dChip analysis and q-values from SAM 
analysis. 
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